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Abstract. Attention has focussed recently on models of inflation that involve a second or 
more fields with a mass near the inflationary Hubble parameter H, as may occur in super- 
symmetric theories if the supersymmetry-breaking scale is not far from H. Quasi-single-field 
(QsF) inflation is a relatively simple family of phenomenological models that serve as a proxy 
for theories with additional fields with masses m ~ 77. Since QsF inflation involves fields in 
addition to the inflaton, the consistency conditions between correlations that arise in single¬ 
clock inflation are not necessarily satished. As a result, correlation functions in the squeezed 
limit may be larger than in single-field inflation. Scalar non-Gaussianities mediated by the 
massive isocurvature held in QsF have been shown to be potentially observable. These are 
especially interesting since they would convey information about the mass of the isocurvature 
held. Here we consider non-Gaussian correlators involving tensor modes and their observa¬ 
tional signatures. A physical correlation between a (long-wavelength) tensor mode and two 
scalar modes (tss), for instance, may give rise to local departures from statistical isotropy 
or, in other words, a non-trivial four-point function. The presence of the tensor mode may 
moreover be inferred geometrically from the shape dependence of the four-point function. 
We compute tss and stt (one soft curvature mode and two hard tensors) bispectra in QsF 
inhation, identifying the conditions necessary for these to “violate” the consistency relations. 
We hnd that while consistency conditions are violated by stt correlations, they are preserved 
by the tss in the minimal QsF model. Our study of primordial correlators which include 
gravitons in seeking imprints of additional helds with masses m ~ 77 during inhation can be 
seen as complementary to the recent “cosmological collider physics” proposal. 
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1 Introduction 

Given that inflation with a single scalar degree of freedom is already able to account for 
observations, one might wonder why we should consider inflationary models with multiple 
fields. One reason is that the flatness required of the inflaton potential calls for a mechanism 
or symmetry to protect it from quantum corrections. This is sometimes described as the eta 
problem. 

A supersymmetric UV completion is often invoked to settle this issue. If supersym¬ 
metry is not broken at energies higher than the inflationary Hubble parameter H, then the 
inflationary vacuum energy will be responsible for breaking supersymmetry, since there is 
no supersymmetric theory in de Sitter space. The theory becomes aware of the difference 
between de Sitter and Minkowski spacetime at a scale that probes effects from the curvature, 
which leads us directly to H. It follows that a natural mass range for fields populating the 
supersymmetric multiplets will not stray too far from m ~ H 

Recently, an interesting proposal was put forward [5] (see also Refs. [6-9]) under the 
name of quasi single-field (QsF) inflation. This paradigm is flexible enough to accommodate 
a number of theories within its parameter space, yet constrained to describe the dynamics of 
a light inflaton field with a number of massive, m ^ H, scalar fields. 

This description sits in a somewhat convenient spot between single-field inflation and 
fully multi-field models, where here “multi-held models” refers to theories with more than one 
light held. Although it is possible that, whatever the mechanism responsible for protecting 
the mass of the inflaton from quantum corrections, it might well do the same for other fields, 
our perspective here will be that of limiting the number of protected fields to the smallest 
necessary, namely one. As we shall see, this choice comes with a number of advantages, chief 
among which stands the more predictive nature of quasi-single held inflation. 

^Again, the underlying assumption here is that supersymmetry is being broken at H. If the breaking 
scale were much higher, then the corresponding massive fields may be integrated out, effectively reducing the 
dynamics to that of a single degree of freedom. Intriguingly, even in the low-energy, single-field regime, one 
may find remnants of a past multi-field dynamics (see e.g.. Refs. [1-4]). 
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A field with a mass of order H decays relatively quickly, so much so that its most promi¬ 
nent (late-time) effects are those mediated by its interactions^ with the light infiaton field. 
It is the “action-only-through-</>” (where (j) is the infiaton) that makes for a predictive theory 
and shields QsF inflation from a host of highly model-dependent characterizations extend¬ 
ing to the post-inflationary evolution. This is in stark contradistinction to what generally 
applies for fully multi-held theories where the extra held(s) directly contributes to late-time 
observables. 

The non-Gaussianities in the scalar sector of QsF inhation have been the subject of 
a thorough investigation in Ref. [5]. Here we will focus on correlations between scalar and 
tensor modes. A correlation between a soft (long-wavelength) tensor mode and two short- 
wavelength curvature huctuations (henceforth tss correlation), for instance, would be instru¬ 
mental in unveiling general features of inhationary models. Even though it may be difficult 
to measure it directly, such a correlation is bound to affect the two-point function of the cos¬ 
mic microwave background (CMB) huctuations and galaxy distributions (i) by introducing 
a quadrupolar correction (if the long-wavelength tensor mode is larger than the size of the 
observable horizon); or (ii) in the form of a local departure from statistical isotropy for sub- 
Hubble long wavelength tensor modes. This local departure from statistical isotropy appears 
in Fourier space as an oh-diagonal correlation between different Fourier modes of the density 
held, or equivalently, a nontrivial four-point scalar correlation function. These features, in 
their interpretation as remnants of held dynamics during inhation, are often dubbed as fossils 
[ 10 - 12 ]. 

A tss correlation signal in standard single-held slow-roll inhation (SFSR), and, more 
generally, in all single-clock models, is partially cancelled by projection (late times) effects, 
leaving an unobservably small (of order /kg) signal [12-14] . This is no longer the case for 
models that violate^ consistency conditions as has been investigated in Ref. [15] for inhation 
with non Bunch-Davies initial conditions, in Ref. [16] for inhation with a non-attractor phase 
[17-19] (see also Ref. [20]), and for solid inhation [21-24]. 

In single-clock models of inhation, consistency conditions (henceforth ccs) [25-28, 30- 
32, 80] relate the n-point correlation functions to (n—l)-point functions in the soft limit of one 
of the momenta. They do not necessarily apply, however, if there are multiple helds during 
inhation. The squeezed limit of cosmological correlators is therefore a powerful instrument 
to probe inhation. 

In this paper we study tss and stt bispectra in QsF inhation in the squeezed limit 
for, respectively, the tensor and the scalar curvature mode. We show that the stt receives 
ccs-violating contributions, thus opening the way to new observables. On the other hand, ob¬ 
servable signatures from tss correlations may not be generated within the minimal realization 
of QsF inhation as introduced in Refs. [5, 33]. We show that, for a violation of ccs involving 
the tss correlation to occur, a non-zero two-point correlation between the soft tensor mode 
and the massive isocurvature mode is necessary. This is forbidden by the symmetries of 
the theory in the minimal QsF scenario. However, such violation may arise in non-minimal 

^Indeed, in QsF inflation, the linear equations of motion for the fields are typically decoupled because the 
would-be mixing term in the Lagrangian is assumed small and can be treated as a perturbation on top of £ 2 . 

^It is important to note here what exactly is to be intended as “violation of ccs”. It is not the non-linearly 
realized symmetry to be violated/broken. Rather, in all that follows, whenever a relation between a “squeezed” 
n-f 1 and an n-point functions will require the knowledge of an additional non-trivial (i.e. not reducible to the 
original two ingredients) n + 1-correlator (typically involving additional fields, see the “extra” in Eq.(4.4)), 
we will say there is a ccs violation. 
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realizations with broken statistical isotropy or in the presence of additional vector or tensor 
modes. 

This work is organized as follows: in Section 2 we review the main features of the model 
and its predictions for the scalar sector; in Section 3 we provide a general discussion on 
signatures from a tensor-scalar-scalar correlator; in Section 4 we present our hndings on the 
three-point correlation functions involving gravitons; in Section 5 we offer comments and 
conclusions. The details of some of our derivations are provided in Appendix A. 


2 The model: review of background and perturbation analysis 


QsF inflation comprises a class of multi-field models characterized by the presence of a light 
held (the curvature mode) driving inhation, and one or more helds (isocurvature modes) 
with masses of order of the Hubble scale H. Curvature and isocurvature modes couple in 
the presence of a turning trajectory in held space [34]. In Ref. [5, 33]'^ a very simple quasi- 
single-held model was analyzed, with one massive held (cr) only, and a minimal coupling to 
gravity. Using polar coordinates in held space, the Lagrangian has the form. 


Sm= j 


\ {d^a){du(T) - Vsr{0) - V{a) 


,( 2 . 1 ) 


where R9 corresponds to the adiabatic direction, parallel to the inhaton held trajectory, and 
a to the perpendicular (isocurvature) direction. For simplicity, the trajectory in held space 
is assumed to have constant radius and to be characterized by a constant angular velocity 
(constant turn). Another simplifying assumption concerns the magnitude of interactions 
between a and 0 at the quadratic level: 

Cc^ = 2a^ReQ5a5e, (2.2) 

where we have replaced R = R + gq, with cjo being the value of a at the minimum of its 
effective potential. 

For small values of the coupling, \9q/H\ <C 1, one can treat Cc 2 as a perturbation 
on top of the second-order Lagrangian for the two decoupled free helds. This allows us to 
separately solve for the linear huctuations and then study the mixing between adiabatic and 
isocurvature modes using perturbation theory. The mixing has been “postponed” up until 
when interactions are considered. The free-held mode functions have the form: 

for m < ^ , (2.4) 

where v = , and m? = V -|- 7 0g is the mass of the isocurvature huctuation^. 

As expected, the larger the mass the faster the mode function of the isocurvature held 

“^See also Refs. [35-38] for further studies on this model and Ref. [39] for the connection between super- 
symmetry and Hubble-mass degrees of freedom (a recent model with e.g. maxion ~ H is found in ]40]). 

®The m > (3/2) H case entails an oscillatory behavior at late times which can suppress the contribution 
from ^CT-mediated correlators. For this reason the authors of [5, 33] focus on the m < (3/2) H case. Note 

however that the lower part of this mass range is still far from the integrating out regime and supports a very 

interesting dynamics [41, 42]. 
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asymptotes to zero after horizon crossing, 


5o-fc(r) ~ , (/cT 0). 


(2.5) 


At tree level, the power spectrum of curvature fluctuations receives a scale-invariant contri¬ 
bution from interactions of the type in Cc 2 (as represented in diagram (b) of Fig.l) that add 
up to the standard power spectrum (diagram (a) of the same figure): 



where C is a function of the mass of the isocurvature field (see Refs. [5, 33] for its numerical 
evaluation) and C ~ —H 56/9q. 


The field a is not in slow-roll (it sits in the minimum of the effective potential) and, as 
a consequence, the magnitude of V can be much larger than the potential in conventional 
SFSR scenarios or in the case of multiple light fields models. As a result, contributions to the 


scalar bispectrum arise from diagrams that are mediated by 5a (as in diagram (d) of Fig. 1) 
that can be large compared to the standard one-vertex diagram (diagram (c) of Fig. 1). One 
can see that the magnitude of diagram (d) would be suppressed if a itself was in slow-roll: 



In the squeezed limit, /cs <C fci ~ A: 2 , and for u ^ Q, the contribution to the scalar bispectrum 


from diagram (d) is found [5, 33] to have the following momentum dependence on the mass 
of the isocurvature fluctuations 



( 2 . 8 ) 


The momentum dependence in Eq. (2.8) is easily understood from the super-horizon decay 


of the isocurvature fluctuations, as given in Eq. (2.4): in models with standard Bunch-Davies 
initial conditions, as is the case for QsE inflation in Ref. [5, 33], a correlation between long 


and short wavelengths can only be generated once the short-wavelength modes approach 


horizon scales; by that time, the amplitude of the long-wavelength mode will have decayed 


by (ri/r3)3/2-- = 


\-u _ 


The shape function in Eq. (2.8) interpolates between local (to which it reduces for 
—)• 3/2) and equilateral (approximately approached for z/ —)• 0). The interpolation between 
local and equilateral shapes (see Ref. [33] for more details) can be easily understood: the 


more massive 5a is, the faster it decays on super-horizon scales and therefore, the larger the 
contribution to non-Gaussianity generated around horizon-crossing scales (equilateral type). 
On the other hand, for lighter masses of 5a, the super-horizon isocurvature fluctuations sur¬ 
vive longer and can contribute to correlations between long and short modes. 

Let us briefly comment on the prospects for detection of the scalar bispectrum in QsF 
inflation. If V /H < 0(1), the leading contribution arises from the interaction in Eq. (2.7). 
The amplitude reads [5]: 



(2.9) 
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(b) 



Figure 1. (a) Standard power spectrum of (. (b) Correction to the power spectrum of ( from 

interactions with a. (c) Bispectrum of C from self-interactions, (d) Bispectrum of C from interactions 
with (7. Dashed line are associated with a, continuous ones with 


While we know that /nl in Eq. (2.9) could be of order 10 (for the largest possible V'” 
values, the constraint being V < H), it is just as true that it could be much smaller (and 
yet dominated by the V contribution). It follows that such a signal may or may not be 
accessible for, e.g., upcoming LSS surveys. 


3 Tensor-scalar-scalar correlation and observational signatures 


A correlation between one long-wavelength tensor and two short-wavelength scalar modes 
manifests itself in the distribution of primordial density fluctuations as an off-diagonal cor¬ 
relation between different Fourier modes of the density field. This is true for any primordial 
field that correlates with curvature fluctuations in the squeezed limit during inflation [10- 
12], but the nature of the correlation depends on the spin of the field. In the presence of a 
long-wavelength mode of the new field, the correlation between two different Fourier modes 
becomes, 


{S{ki)5{k2)),^p^g^ = (27r) 


h)h;(KHl(K)kuhi 


i:(3) 


(3.1) 


(3) 

where M' is the Fourier transform of the field, p its polarization, <5- ' - indicates the delta 

. ..kn 

function of argument ki + ... + kn- The function fp is related to the three-point correlation 
of the field with the curvature fluctuations by 

Bp{K, kiM) = Pp{K)fp{kiM)e\.{k)kukj = BpiK, ki, k2)e%{k)kuk2j, (3.2) 


where Pp is the power spectrum. Global statistical isotropy requires that / is only a function 
of the norm of ki and k 2 and of the angle between the two vectors. 
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Figure 2. (e) Representation of a graviton-exchange diagram in the four-point function (which we 
mimic in the soft K <^ki limit). Note the black shaded area stands for a generic type of interaction, 
(f) A pictorial representation of the momenta configuration of the non-trivial four-point function we 
are effectively probing in Eq. (3.4). 


Whenever the long-wavelength mode is super-Hubble, ki and ^2 are indistinguishable 
from one another. In this case, the two contributions on the right-hand side of Eq. (3.1) can 
be condensed into a single diagonal term that corresponds to the sum of the regular scalar 
power spectrum plus an anisotropic correction. 


P{k) ~ P{k) 


1 Pa'yijkikj 


(3.3) 


where is the tensor perturbation in our observable volume. Here a is related to the 
amplitude of the tss correlation, a ~ B{kL, ks, ks)/P{ks)P-y{kL). Eq. (3.3) shows that the 


expectation value { S{k) ) about the direction K of the long-wavelength mode can be 


isotropic (if the extra field is a scalar) or anisotropic if the extra field is a vector or tensor 
field [10, 11]. In this work, the field is the tensor mode from the metric ( 7 ^). 

Null searches for the quadrupolar anisotropies described by Eq. (3.3), when the new 
field is a tensor field, both in the CMB [43-48] and LSS [49, 50], have resulted in a lower 
bound < 0 ( 0 . 01 ) on the magnitude of the quadrupole. 

If the long-wavelength mode is sub-Hubble, any set of two scalar modes of wavenumbers 
ki and k 2 with ki + k 2 + K = 0 can be used to estimate the amplitude of the tensor mode, 
providing a powerful probe that may well complement other existing ones aimed at the search 
for primordial gravitational waves (including measurements of the B-mode CMB polarization 
[51-56], gravitational lensing effects in the CMB [37, 57-59], LSS [59-63] and 21cm cosmology 
[64-66], and direct gravitational-wave searches [67-74]). Eor a stochastic GW background, 
the minimum-variance estimator for the tensor amplitude is [ 10 ], 








(3.4) 


where P/ = P^{K)/A^ is a fiducial power spectrum for the tensor modes. The quantity 
%{K) is the optimal estimator for the amplitude of a single Eourier mode. 


%{K)^P^{K)^ 

k 


2yptot(^)ptot(|^_^|) 


(3.5) 
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and Pp{K) its variance 


Pp{K) = 


\Bp{K,k,K-k)/P^{K)\^ 
2^ 2VP^^\k)P^^WK-k\) 


(3.6) 


In the above formulas V = (27r/A:niin)^ stands for the volume of the survey and is the 
total scalar power spectrum, including signal and noise. The variance of the estimator in 
Eq. (3.4) 


^7 ' = 2 5 W) > (3-7) 

K,p 

may be used to estimate the smallest GW amplitude that can be detected for a given survey 
size. 

As already mentioned, the actual observed quantities are found after an additional step. 
Specihcally, one needs to subtract from the primordial tss correlation a term accounting for 
late-time effects [12-14], 

Sobsiki, ks, ks) = B{kL, ks, ks) - BcdkL, ks, ks), (3.8) 

where the tss correlator has the generic form, 

k2, ks), (3.9) 

and we set /cl = fei <C k 2 ,ks ~ ks- The second term on the right-hand side of Eq. (3.8) 
corresponds precisely to the expression, to leading order in powers of (ks/ks), for the squeezed 
limit tss correlator as dictated by consistency conditions [25]; 

Bcc{kL,ks,ks) = ■ ( 3 - 10 ) 

In models where the consistency conditions are satisfied, the leading-order part of the pri¬ 
mordial signal would then be, for the most part, cancelled by projection effects, leaving an 
unobservably small 0{k‘j^/kg) signal. As a result, observations (either direct or indirect) of 
a tensor-scalar-scalar correlation would rule out single-clock models of inflation. 

Eq. (3.4) and (3.5) reveal that we are effectively probing a specialized four-point function 
in the limit that approaches the so-called counter-collinear conhguration, see Fig. 2(e). 

4 Correlators with tensors in QsF inflation 

Similarly to the scalar bispectrum, the tensor-scalar-scalar three-point function may be gen¬ 
erated by direct interaction of tensor modes with scalar curvature (as in diagram (g) of 
Fig. 3) 


BC5 


aR^e^ 

2~hP 


^ijdiCdjC ) 


(4.1) 
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(g) 



(h) 




Figure 3. (g) Tensor-scalar-scalar correlation from usual coupling between 7 and C. (h,i) Leading 
order tensor-scalar-scalar correlator mediated by a. Dashed line are associated with a, wiggly lines 
with 7 and solid lines with 


or as the result of interactions between the external fields and the isocurvature mode (as in 
diagram (h) and (i) of Fig. 3). Diagram (h) arises from the interactions between 6a and the 
graviton, as well as from the interactions between 6a and 

Cc, = -2 C, Cc, = -^ijdi6adj6a . (4.2) 

The diagram can be computed using the in-in formalism [75]. We refer the reader to Ap¬ 
pendix A for the details of the derivation and report here the final result in the squeezed 
limit 






42 / 
^0 / 


= w{u)^Py{kL)P^^Hks) 


'm 


(4.3) 


Here, P'y{k) = {H‘^)/{Mpk^) is the tensor power spectrum, and P^^\k) = {H‘^){2P?6Qk^) 
the leading order part of the scalar power spectrum. In order to maintain analytical control, 
we have calculated the tss correlator directly in the squeezed limit. One could schematically 
think of the complete tss correlator as, for example, the sum of two contributions, one that 
dominates in the squeezed limit and the other that dominates, e.g., in the equilateral limit. 
However, our variance in Eq. (3.4) accounts for a hierarchy between momenta associated 
with tensors and scalar modes, it is built from an observable evaluated in the squeezed 
limit, therefore one would expect that it would not be severely affected by contributions 
from non-squeezed momenta configurations. The function w{v) is shown in Fig. 4 for a 
range of v values. We defined ki and ks as, respectively, the (tensor) long- and (scalar) 
short-wavelength modes. 

The tss correlation here is of the local type: no suppression as that appearing in Eq. (2.8) 
is now present because 6a does not effectively carry the soft momentum. One can see this 
at the level of the Feynman diagrams for the IN-IN formalism employed in our calculations. 
Whilst in Fig. 1(d) one isocurvature fluctuation 6a is necessarily contracted with an external 
scalar fluctuation, forcing its wavenumber to be soft due to momentum conservation, this is 
not the case of Fig. 3(h), where the external tensor attaches directly to a three-point vertex. 
Correlations in the squeezed limit are generated as soon as the isocurvature modes have 
reached horizon size, before they undergo any damping. 
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Figure 4. Plot of numerical values of the coefficient w(i/) introduced in Eq. (4.3) for i/ ranging from 
0.06 up to 0.4. 


Notice that the r.h.s. of Eq. (3.10) can be expanded to next-to-leading order in slow- 
roll: in this case it will include contributions proportional to {9/H)‘^ (see e.g. the correction 
from the isocurvature mode to the curvature power spectrum in Eq. (2.6)). This observation, 
along with our findings in Eq. (4.3), are suggestive of the fact that consistency conditions 
may be preserved® for the tss correlator and they should be carefully checked. 

This can be done by computing diagram (i) of Fig. 3. It is easy to verify that, modulo 
an overall numerical coefficient (in the form of some function, g{i 2 ), of the isocurvature field 
mass), diagrams (h) and (i) have the same parametric and momentum dependence. One 
may then follow the steps outlined in Appendix A to check the consistency conditions with 
a straightforward but lengthy calculation. Given the approximations and the final numerical 
integration involved in our squeezed-limit computation, this may not be the most convenient 
route towards verifying an exact cancellation. We opt instead for a more general procedure 
that allows us to also identify what is to be expected of a given model for the consistency 
conditions to be violated. We outline this in the following subsection. 


4.1 Consistency conditions for inflation with extra flelds 

We will be following the notation and logic of [78-80] to show that the general form of the 
tensor-scalar-scalar correlator, to leading order in the squeezed limit, is as follows 


lim - — 

fci^o -^7(^1) 



3 


a=2 


d 

dki 


{XX)c + Mra], 


(4.4) 


where a contribution proportional to the first one on the right-hand side is the only one found 
when, in our nomenclature, the model preserves consistency conditions '. We will briefly dis¬ 
cuss the origin and form of the extra contribution that appears when these are violated. 

®We are grateful to L. Bordin, P. Creminelli, J. Norena for illuminating discussions on this point. 

^Here [extra] stands for everything that makes up the relation dictated by the symmetry. It is important 
to stress that, in general, only a subset of [extra] violates the ccs. The culprits are in particular quantities, 
such as some of those contributing to the last term in Eq. (4.13), which cannot be trivially reduced to the 
usual (l.h.s. and r.h.s.) two terms in Eq. (4.4) above. 
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In Eq. (4.4), a prime stands for an (omitted) momentum-conserving delta function and a 
subscript “c” indicates connected diagrams. We will drop the connected diagram subscript 
from here on. 

First off, one needs to identify the (non-linearly realized) symmetry that generates the soft- 
tensor consistency conditions. This corresponds to an anisotropic spatial rescaling of coordi¬ 
nates 


X* —)> X* -|- Sij , 


(4.5) 


where is a symmetric, transverse and traceless tensor. Under this coordinate transforma¬ 
tion the curvature fluctuation shifts linearly (<5a;C ~ Sijx^diC), whereas crucially the tensor 
modes transformation has a non-linear contribution {5u}'yij ~ Sij). Introducing the Noether 
charge, Quj, associated with the symmetry transformation, one hnds: 

{Sco'yiji^) = i ( [Quj,7ij{^] ) ~ (27r)3(l(3)(^5y . (4.6) 


Consider now the effect of the transformation on the operator while, at the same time, 

introducing a complete set of states : 


n,k 


n,k 


(4.7) 


We define one-particle states generated by C, 7 and by the orthogonal projection of the 
isocurvature mode a as: 


17 




|i P 






(4.8) 


where 




C? - 


-1 




7; 


ij,k 


(4.9) 


and we have used the fact that (C 7 ) = 0 by rotational invariance. Using this in Eq. (4.7) we 
obtain: 


CkS, 


k2^kz 


z ^Im Cfcj)(CfcjCfc2Cfc3) + P-f ^{^^)iQ^lij,kx)^'^ij,kSk2^i 


ks' 


ki 




(4.10) 


The l.h.s. of the equation just above may be obtained by simply writing explicitly the linear 
part of the transformation for This then gives precisely the well-known leading order 
contribution [25] to the tensor-scalar-scalar correlator (hrst term on the r.h.s. of Eq. (4.4)) 


( Q^^^kz^kz ) 



(4.11) 


Moving on to the r.h.s. of (4.10), the first contribution, proportional to {QuC^^ vanishes 
because C, transforms linearly under the tensor symmetry. Using Eq. (4.6), the second con¬ 
tribution gives 


fci 


Im 


P^ ~' Ji™ ^7 \^^)Si3{iij,kSk2kz) 


(4.12) 
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which corresponds to the l.h.s. of Eq. (4.4). A violation of the consistency conditions can 
therefore only arise from the third contribution to Eq. (4.10). This is given by 


)(dr Cr Cr;J ~ (Qa;^) - 


{QojOiCcr) {Qujl){l(y) 


Pc 


P. 




(4.13) 


Notice that the dots in Eq.(4.13) stem from our expanding the state |(t) as a function 
of the more “physical” |ct); the expansion is regulated by the small parameter e(z^)^ = 
\{C,a)\‘^/ {P(^ Pa) ■ The first two terms on the r.h.s. of (4.13) are again equal to zero from 
the linear action of the tensor symmetry on Q and a. The third contribution is zero in the 
case of QsF inflation because rotational invariance requires (uy) = 0 . This goes to show 
that a ccs-violating contribution to the tss correlation can be produced from a non-zero 
tensor-isocurvature mode correlation. While it is apparent that Lorentz invariance forbids 
any direct tensor-isocurvature quadratic coupling to arise from the Lagrangian in (2.1), the 
validity of our argument also applies at loop level: the absence of any source breaking ro¬ 
tational invariance in the model ensures that the tensor-isocurvature correlation is equal to 
zero. Symmetry considerations indeed require the following 

^ ^^^Hk + p)P{k)Zij , (4.14) 


where Zij is a tensor, independent from the metric and the wave vector k, enjoying the same 
properties as 'jij (symmetric, traceless, transverse). A consequence of this result for QsF is 
that also diagrams as the one in Fig. (5), which one would naively expect to be only mildly 
suppressed w.r.t. tree-level contributions®, are null for symmetry reasons. 

To conclude, a detection of a tss correlation in the limit of a soft tensor mode may 
point towards massive isocurvature fields during inflation, but not in the minimal set-up 
analysed in this work, where the massive isocurvature fields are scalars and the background 
is isotropic. One may break statistical isotropy by introducing vector degrees of freedom. 
It would also be interesting to explore scenarios with additional tensor modes. It is worth 
investigating possible extensions, along these lines, of the minimal set-up of QsF inflation, 
and assess whether the ensemble of constraints on the theory would allow for an observable 
fossil signal. 


Within the minimal QsF model, it is also interesting to ask whether other correlators 
involving tensor modes may produce a violation of ccs. Consider for instance the scalar- 
tensor-tensor (stt) bispectrum. In the limit of a soft curvature fluctuation, one can follow 
the arguments outlined in this section for the tss correlation. These show that a necessary 
conditions for the stt correlation to violate ccs is that of C being sourced by a long-wavelength 
isocurvature mode. This is the case for the diagrams represented in Fig. 6. A typical 

^Before relying on symmetry arguments, let us consider the case when a self-couplings are not much smaller 
than unity, that ist^ /H < 0{1). One may expect diagrams with isocurvature fluctuations running in the 
loop not to be significantly suppressed in comparison with tree-level contributions such as those in Fig. 3(h,i). 
The diagram in Fig. 5, for instance, would naively have a squeezed limit amplitude 





P.,(fci)Pc(fe2) 



(4.15) 


However, for such a diagram to be non-zero one would need to introduce some degree of anisotropy in the 
model. 
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Figure 5. Example of loop diagram contribution to the tss bispectrum. In QsF inflation any diagram 
that invoves a correlation between one tensor mode and one isocurvature mode is equal to zero for 
symmetry reasons. 


contribution to the latter will display a {ki/ks)^^ dependence in the squeezed limit which 
is not captured by the standard ccs. 

The symmetry of interest for the squeezed stt correlator is a dilatation x* —)> (1 + A) x*, 
under which the curvature fluctuation transforms asC C + -^(l + ^‘ dxO- From the 
transformation of tensor modes under dilatation, 5d 'Jij ~ A , one finds: 




) = ( 


'k2 'ks 


i A 


3 + 'y ^ ka 

a=2,3 



'k2 'k3 


(4.16) 


which corresponds to the usual expression for the squeezed stt correlator [25]. Similarly to 
Eq. (4.7), one can expand the l.h.s. of Eq. (4.16) and introduce a complete set of states 


( (< 3 rfKfc)(^fcl 7 fc, 7 fc 3 > - {VkJ^k){^k\Qd) 

n,k n,k 

As before, we consider one-particle states for C, 7 and a. One hnds 


[Qq,7^ 


i Im 


P^^Qd 0 (C Y) + Py'iQd 7) (7") + Pd'iQd d) (d 7 ") 


(4.17) 


(4.18) 


The (QdC) term gives a contribution proportional to iQdj) term is zero 

from tensor modes transforming linearly under dilatation. A non-zero contribution originates 
from the {QdY term, in the presence of a correlation between isocurvature and curvature 
modes, which is the case for QsE inflation. Schematically this term gives 


{QdYYY) ~ [{Qdcr) - ^{QdOiCY - P^ ^{QdYi'yY) YY) 

~ \i5{^P^^{CYYY) H— , (4-19) 

where the first term in the first line on the r.h.s. of Eq. (4.19) is again zero because a itself, 
as a scalar and in general a non-gravitational degree of freedom, transforms linearly under 
the symmetry. The final result is that the consistency-condition violating contributions for 
the squeezed stt originates from: 


P, 


-1 


(%d£7g) ^ {y,YIi^ 


^ 2^73 
k2 


(4.20) 


and in Fig. 6 we represent some of such diagrams. 

It would be interesting to identify possible signatures of stt correlation functions in QsF 
inflation and in more general models with fields with m ~ 0{H). We leave this for future 
work. 
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Figure 6. Representation of ccs-violating contributions to the stt correlation. The soft scalar curva¬ 
ture mode is directly sourced by a (soft) isocurvature mode. 


5 Conclusions 

The combination of CMB and LSS data at our current (or near-future) disposal represents 
an unprecedented source of information on the physics of the early Universe. Most of the 
cosmological evolution probed takes place at very high energies, much higher than those 
one can hope to reach through particle colliders. Cosmological probes then might well be 
our best chance at an observational window on UV phenomena. Naturally, unravelling the 
cosmological information and, most notably, mapping it to operators in the Lagrangian of 
the underlying theory, is a complicated endeavour. 

Many successful efforts in this direction already populate the literature, and Ref. [42] 
recently systematically classified the properties of inflationary curvature correlation functions 
in the presence of additional fields with Hubble-scale masses (see also [77] for a very general 
analysis on inflationary correlators in the presence of massive fields). The oscillatory (or, 
alternatively, power-law) behavior of correlators in the squeezed limit provides a measure of 
the mass of the fields while the angular dependence encodes knowledge about their spin [10] . 
Here as elsewhere the study of the squeezed limit of correlations functions proves crucial 
for investigating new degrees of freedom during inflation. One may well see part of the 
present work as complementing the one in [42] in that we study the effect of massive fields 
on correlators that also contain gravitons. 

Here we investigated the observability of squeezed three-point correlation functions in¬ 
volving both curvature and tensor fluctuations in QsF inflation. The latter is a phenomeno¬ 
logical inflationary scenario comprising one or more fields with mass of order Hubble. In 
order to assess the observability of such signals, the first step is to verify whether these bis¬ 
pectra violate single-clock consistency conditions. We checked from general principles that, 
in QsF inflation, this is the case for correlations among one (soft) curvature mode and two 
tensor modes (stt). The latter opens up the possibility of new late-time observable effects in 
CMB fluctuations or in the galaxy distribution. It would be interesting to further investigate 
such a possibility in this as well as other inflationary models that include m ~ 77 fields. 

An observable of great interest, as pointed out in previous works in the context of other 
inflationary models [16], is the tss correlation (one long-wavelength tensor and two scalar 
modes). In [10]-[12] it was shown how, in the presence of such a primordial correlation, a soft 
super-Hubble tensor mode would induce a power quadrupole, whereas a sub-Hubble tensor 
mode would introduce local departures from statistical isotropy in the form of an off-diagonal 


- 13 








correlation between different Fourier modes of the density field. The latter could be used as 
a probe for the amplitude of tensor modes and, for a given value of the Hubble parameter 
during inflation, one can estimate the size of the survey necessary to access the physics at that 
scale H. This is not the most direct way of probing the existence of tensor modes. Rather, 
it is a specific probe that relies on, and therefore crucially contains information about, the 
existence of a non-Gaussian signature arising from the multi-field and non-single clock nature 
of the QsF model. It is in other words a route to seeking critical information on the nature 
of the inflationary dynamics. 

In this paper, we show that violation of consistency conditions for the tss correlation 
would occur in the presence of a non-zero two-point correlation between the tensor mode 
and the massive isocurvature. Symmetry requirements (statistical isotropy) in the minimal 
QsF inflation realization considered in this paper forbid such a coupling. However, tensor- 
isocurvature correlations may well be generated within non-minimal extensions of QsF. This 
may be the case if one breaks rotational invariance e.g. by extending the particle content 
to include additional vector or tensor degrees of freedom. These are interesting possibilities, 
worth exploring. 

One may also wonder how these results may be extended to the fully multi-field case 
(in the sense of non-decaying light degrees of freedom alongside the inflaton field). As men¬ 
tioned, a high degree of model dependence is intrinsic to these setups. Nevertheless, the 
main properties of QsF inflation we have employed here, also characterizes the dynamics of 
multiple light fields. We leave this to future work. 
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A Computation of ^cr-mediated tss correlation in the squeezed 
limit 


The expectation value of an operator 0 at time t can be computed using the in-in formula, 


(0(Q)= 0 


T exp i dt Hj{t ) 


'to 


0/(t) 


" y 'f s 


T exp '—i dt Hj{t ) 


fto 


where Hj is the interaction Hamiltonian. For a three-vertex diagram one finds 

rt rti rt2 


(0(Q) = -2X 


dti f dt2 f dts {Qi{t)Hj{ti)Hi{t2)Hi{t3)) 

to J to J to 

t 


+ 21 


dii f dti f dt2{Hi{ii)&i{t)Hi{ti)Hj{t2)) 

Jto Jto 


0) , (A.l) 


(A.2) 


(A.3) 
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where X stands for the imaginary part. 

For an operator 0 ~ and an interaction Hamiltonian Hj = H 2 + H 3 , where 

H 2 ~ Sa59' and H 3 jijdidadjda, from the right-hand side of Eq. (A.2), three terms arise 
(the same goes for Eq. (A.3)). The number of permutations for each of the six different terms 
is 4, so there is a total of 24 terms. The right hand side of Eq. (A.2) is then given (modulo 
an overall factor) by the sum of the following terms: 


-2eijk2ihj 'yk^{0)uk2{0)uk:^{0) x X 


f 

' —00 
/*X2 


dTia^{Ti)-il^{Ti)vk 2 {ri)vk 3 {Ti) (A.4) 


X / dT2a;\T2)vl {T2)u^AT2) / (ts) uF (rs)^^* (ts) 


' —00 


-|- 3 perms. 


-2e,^ A:2i/c3j7fci(0)Mfc2(0)«fc3(0) x 


dria^(ri)u. (ri)ufc (n) 


rri 


rr2 


X drao {T 2 ) jl^{T2)v*k^{T2)vkAT2) / dna (ra) (rs)^^ (ra) 


(A.5) 

-|- 3 perms. 


-2e-)X2ifcaj7fci(0)'Ufc2(0)ttfc3(0) x X 


dna^{Tl)Vk^{Tl)Uk^{Tl) 


X / dT2a^{T2)vk^{T2)uk^{T2) / {n) 


Eq. (A.3) is given by the sum of the following terms 


(A.6) 

-|- 3 perms. 


2e^jk2ik‘ij-il^A))uk^A))uk^{A) x X 


dfi a^(ri)7fci(H)^'fc2(H)%(H) 


X / dTia^{Tl)vl^{Tl)u^^{Ti) [ dT2a^{T2)v*k^{T2)u^^{T2) 

J —00 J —00 


-|- 3 perms. 


2efjk2ik3j Jki {0)ul^ {0)uk^ (0) x X 


(ifia^(fi)u. (fi)m (fi) 


(A.8) 


X f dTia‘^{Ti)jl^{Tl)vl^{Tl)Vk^{Tl) f dT2a^{T2)v*k^{T2)Uk^{T2) 
J —00 J —00 


+ 3 perms. 


2e^jk2ik3j 7 fc^ (O)^^^ (0)?Xfc3 (0) x X 


dfia^{fi)vk^ifi)uk2{fi) 


X / dria^(ri)ufc3(TiX3(ri) / dT2 {T2) Ak^{T2)vl^{T2)vl^{T2) 


(A.9) 

-|- 3 perms. 
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In the equations above one uses the Fourier-mode decomposition, 

Sk 


(fk 


-ik-x 


atUkir) + a} 


-k 


Sa{T,x) = j "^'^\bkVk{T) + bKvl{T) 




(27r)^ 

d?k 

(2vr)' 


-ik-x 


E4(b 


A=± 


4lkiT)+ ic^r) -fk{x) 


for the fields. The mode functions are given by 


Ukir) ={I + ikr) e , 7k{r) = ^ 


R^/W 


Mp\/W’ 


(1 -|- ikr) e 


—ikr 


Vk 


(r) = Hi^\-kT), 


(A.IO) 

(A.ll) 

(A. 12 ) 

(A.13) 

(A.14) 


where i' = — {m/HY and m?/H"^ <9/4. Consider the first permutation of Eq. (A.4) 

and replace the expressions for the mode functions. One then finds, 


- etkoik 


7r2 ^2 1 


b 2i 3j 25 klk2k3 


xl 


fO l-Tl |•T2 

dn j dT2 / drs 

J —OO 


i-n) 


f —OO J —OO 


(-T2)^/2(_7-3)1/2 


(A.15) 


(1 - ffein) 


Working in the squeezed limit ki k 2 , k^, it is convenient to perform a change of variables 
introducing Xi = ksn. Then the result in Eq. (A.15) becomes 


- ehkoik 


TT 




1 


ij 2i 3i 25 ^2 ^4 klk2kl, 


xX 


rO rxi pX 2 

/ dxi / dx2 / dx^ 

'—OO J—OO J—OO 


{-Xl) 


(-X2)1/2(-X3)V2 


(A.16) 






Xl 




X2 


e*"3i/(2) [_^3] 


Next, one can use the approximation k 2 — k^^ (valid in the squeezed limit) in the first and in 
the third Hankel functions of the previous expression as well as in the second exponential. 
The following approximation is also allowed 




~ 1 . 


(A.17) 


To see why notice that ki ^ 3 . Therefore in order to have (A:i/A: 3 )|xi| > 1 one would need 
|xi| ^ 1, but since xi is the argument of Hankel functions in Eq. (A.16), then in the limit 
of very large |xi| these functions would be rapidly oscillating and suppress the value of the 
integrals. 

With these approximations, the final result for Eq. (A.16) then becomes 


- 4jk2ik3j 


7r2 


1 


25 M|,i?4 A:3A;3 


i rO rxi r-X2 

/ dxi / dx 2 / dx 3 

J—OO J—OO J—OO 


i-Xl) 


(-X2)^/2(_a;3)l/2 


(A.18) 


e“3/^(2) [_a;3] 
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where unit vectors k = k/k have been introduced. With a similar procedure, the first 
permutation in Eq. (A.7) gives 


7 r 2 772 I 


eijk2ik^j 25 jv 72^4 


X X< 


/ U p\j pxi 

dxi / dxi / dx2 

-00 J —00 J —00 


i-xi) 


-00 (-Xi)V2(_^2)1/2 

2 


(A.19) 




The total {'y6656) is obtained by adding up all the permutations listed in Eqs. (A.4)-(A.9). 
Unlike what happens for the {596959) correlation in the squeezed limit [5], for the tss cor¬ 
relator one expects that all of the permutations will have the same momentum dependence 
as in Eqs. (A. 18) and (A.19); e.g. ~ {ki and ks being respectively long and short- 

wavelength modes). 

Summing up all permutations in Eqs. (A.4)-(A.6), one then finds, 


^inksiks 


^2 

MlR^klkl 


xXi 


dx^ 


/ Xi |■X2 

dX2 / 

-00 J —00 


i-xi) 


(-X2)V2(-X3)V2 

i-X2) 


(A.20) 


X 1 -^ 2 ] [-13] + 

xfllD [-x,] + [-n] (iJ<"'|-i'3l)" 


The sum of Eqs. (A.7), (A.8) and (A.9) give 


-^ijksiks 


7r2 772 


1 


^ 23 m 2 7?4 klkl 


xx| / dxi j dxi 
I J —00 J —00 


rxi 


dX2 


i-Xl) 


(-Xi)1/2(_x2)V2 


-Xl 


IXi 


[-ij e“3H(2) [-12] + ■ , [-ill Hl'’ [-n] [- 12 ] 

(-Xi)T4(_x2)V4 


+ 


[-ill [-Ill (ifi='l-i2])' 


(-Xl 


(A.21) 
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